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NOMENCLATURE 


matrix of coefficients of equations of motion 

control matrix in equations of motion 

gust input matrix in equations of motion 

system transfer function 

quadratic performance index 

linear state variable feedback gain matrix 

gust correlation length 

wing torsion degree of freedom 

wing vertical bending degree of freedom 

wing chordwise bending degree of freedom 

matrix of weights on state variables in quadratic performance index 
(normally diagonal) 

matrix of weights on control inputs in quadratic performance index 
(normally diagonal) 
rotor radius 
gust spectrum 

spectrum of state variable response to gust, Sg • |h| 2 S w 

gust input vector 

control vector 

forward speed 

state variable vector 

blade flap degree of freedom; or amplitude of tip path plane tilt 
response, C3j c 2 ♦ Bis 2 ) 1 ^ 2 
rotor coning degree of freedom 
rotor tip path plane pitch degree of freedom 
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B ls rotor tip path plane yaw degree of freedom 
y blade Lock number 

6p wing flap input 

63 pitch/flap coupling 

C blade lag degree of freedom; or amplitude of cyclic lag response, 

(C lc 2 * O l/2 

C damping ratio of eigenvalue, - ReX/ J X | 

C Q rotor collective .lag degree of freedom 

i 0 time derivative of c 0 — rotor speed perturbation for windmilling rotor 
case 

C lc rotor cyclic lag degree of freedom 

C ls rotor cyclic lag degree of freedom 

© 0 rotor collective pitch input 

0 lc rotor lateral cyclic pitch input 

0j s rotor longitudinal cyclic pitch input 
A eigenvalue (root) of system 

Vg rotating natural frequency of blade flap motion 
rotating natural frequency of blade lag motion 
o w rms gust velocity 

Og rms response of state variable 

0| rms velocity of state variable 

o“ rms acceleration of state variable 

(1) frequency (in transfer function or spectrum) 

ft rotor rotational speed 

(*) time derivative 
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SUMMARY 

Optimal control theory is used to design linear state variable feedback 
to improve the dynamic characteristics of a rotor and cantilever wing repre- 
senting the tilting prbprotor aircraft in cruise flight. The response to a 
vertical gust and system damping are used as criteria for the open and closed 
loop performance. The improvement in the dynamic characteristics achievable 
is examined for a gimballed rotor and for a hingeless rotor design. Several 
features of the design process are examined, including: using only the wing 

or only the rotor dynamics in the control system design; the use of a wing 
flap as well as the rotor controls for inputs; and the performance of the sys- 
tem designed for one velocity at other forward speeds. 


INTRODUCTION 

The tilting proprotor aircraft is a promising concept for short haul, 
V/STOL missions. The successful application of this concept will require an 
aircraft with good ride qualities. The combination of large flapping rotors 
operating at high inflow ratio on the tips of flexible wings leads to 
increased response to atmospheric turbulence (compared to a conventional air- 
craft of equivalent size), and hence increases the desirability for an auto- 
matic control system to improve the dynamic characteristics. This report is 
an examination of the basic features of the control systems and design tech- 
niques required for this vehicle. 

To represent the tilting proprotor aircraft in cruise flight, a mathe- 
matical model for a rotor in axial flight on the tip of a cantilever wing is 
used. While the design of actual control systems will of course require a 
model of the complete aircraft dynamics, the proprotor and cantilever wing 
model includes the basic features of the dynamic behavior and hence is satis- 
factory for an investigation of the basic control system characteristics. 
Optimal control theory is used to design linear state variable feedback to 
minimize a quadratic performance index. As a criterion for the open and closed 
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loop performance of the system, the response to a vertical gust is examined. 
Of specific interest are the flap and lag response of the rotor blades, and 
the vertical acceleration of the wing. The damping of the wing modes is also 
a critical indicator of the dyne tic characteristics. The improvement of the 
dynamic characteristics achievable with such control systems is examined for 
two proprotor designs — a gimballed, stiff- inplane rotor and a hingeless, 
soft-inplane rotor. Then several features of the design process are investi- 
gated, including: the effect of considering just the rotor or just the wing 

dynamics in the control system design; the use of a wing flap as well as the 
rotor controls for inputs; and the performance of the system designed for one 
velocity at other cruise speeds. 


PROPROTOR AND CANTILEVER WING THEORETICAL MODEL 

The theory developed in reference 1 will be used for the mathematical 
model of the proprotor and cantilever wing dynamics. 

Figure 1 shows the proprotor configuration considered for the theory. 

The rotor is operating in high inflow axial flight on a cantilever wing. This 
configuration incorporates the features of greatest importance to the aircraft: 
the high inflow aerodynamics of a flapping rotor in axial flow and the coupled 
dynamics of the rotor/pylon/wing aeroelastic system. Many features of the 
aircraft-coupled wing and rotor motion may be studied with such a model, 
theoretically and experimentally, with the understanding that the model must 
eventually incorporate the entire aircraft. 

The theoretical model of the proprotor consists of nine degrees of free- 
dom: the first mode flap (out of disk plane) and lag (inplane) motion for 

each of three blades; and vertical bending, chordwise bending, and torsion for 
the cantilever wing. The degrees of freedom of the individual rotor blades 
are combined into degrees of freedom representing the motion of the rotor as 
a whole in the nonrotating frame. Thus the rotor flap motion is represented 
by the tip path plane pitch and yaw (8 lc and 6j s ) and coning (8 0 ) degrees of 
freedom. The rotor lag motion is represented by cyclic lag (c lc and c, s , the 
lateral and vertical shift of the rotor net Cw.»cer of gravity) and collective 
lag (; Q ). Wing vertical and chordwise bending (qj and q 2 ) and torsion about 
the elastic axis (p) complete the nine degrees of freedom. 

The rotor blade motion is represented by first mode flap and lag motion 
and is assumed to be pure out-of-plane and pure inplane deflections, respec- 
tively, of the blade spar. For the gimballed and hingeless rotor blades con- 
sidered here (except for the flap mode of the gimba-led rotor), there is, in 
fact, some elastic coupling of the flap and lag modes, so that both out-of- 
plane and inplane motion participate in each mode. In the coefficients giving 
the aerodynamic forces on the rotor, it is further assumed that the mode 
shapes are proportional to the radial distance from the hub, i.e., equivalent 
to rigid body rotation about a central hinge. The model based on these two 
assumptions, which considerably simplify the aerodynamic and structural terms 
of the rotor equations, proves to be an adequate representation of the funda- 
mental proprotor dynamics. 
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The theoretical results presented here will be for the rotor operating 
unpowered, i.e., windmilling or autorotation operation. An important element 
of autorotation dynamic behavior is the rotor speed perturbation. With no 
restraint on the rotor shaft rotation, this degree of freedom has considerable 
influence on the aeroelastic behavior of the proprotor and wing. The rotor 
speed perturbation is modeled by using the collective lag mode c 0 . By set- 
ting the rotating natural frequency of this mode to zero, i.e., no spring 
restraint, 5 0 becomes equivalent to the rotor speed perturbation (the natural 
frequencies of the cyclic lag modes, 5 lc and c ls , are not set to zero). 

The proprotor operating in high inflow has simpler aerodynamics than the 
helicopter rotor in forward flight. As in the case of low inflow (i.e., the 
hovering helicopter rotor) , the symmetry of axial flow results in a corre- 
sponding symmetry in the equations of motion; it also means that the equations 
of motion have constant- coefficients. In high inflow there is the additional 
fact that both out-of-plane and inplane motions of the blade produce signifi- 
cant angle-of-attack changes at the sections, and the resulting lift increment 
has significant components both normal to and in the disk plane. Hence the 
rotor aerodynamic forces are primarily due to the lift changes produced by 
angle of attack changes, i.e., the C 5 terms in the aerodynamic coefficients. 

Gt 

This is in contrast to low inflow, where, for example, the inplane blade 
motion produces significant contributions to the forces by the lift and drag 
increments due to the dynamic pressure changes, i.e., the and Cj terms 
in the aerodynamic coefficients. As a result, high inflow aerodynamics are 
well represented by considering only the C£ forces. If, in addition, the 

lift curve slope is assumed constant, then tfie aerodynamic coefficients depend 
only on two parameters, the Lock number y and the inflow ratio V/ftR. 

This nine degree of freedom model will have nine roots or eigenvalues 
(really nine pairs of complex roots) and nine corresponding eigenvectors or 
modes. Each mode involves motion of all nine degrees of freedom. The modes 
are identifiable by their frequencies (which are near the nonrotating 
uncoupled natural frequencies for the rotor modes) , and also by the partici- 
pation of the degrees of freedom in the eigenvector. The nine modes will be 
denoted as follows (the approximate uncoupled, nonrotating frequency of the 
mode is given in parentheses) ; 

p wing torsion (top) 

q 2 wing vertical bending (“q^ 

q 2 wing chordwise bending (“q 2 ) 

6 coning (vg) 

S+l high frequency flap (vg + n) 

6-1 low frequency flap (vg - £1) 

5 collective lag (v ? ) 

5+1 high frequency lag (v^ ♦ $1) 

5-1 low frequency lag (v;. - n) 

The basic theoretical model will consist of all degrees of freedom, 
autorotation operation, and just the rotor aerodynamic forces. The wing 

aerodynamic forces are also included, based on a strip theory calculation. 
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The theory described above will be applied to two full-sea^ proprotors. 

The first is a 7.6 m (25 ft) diameter gimballed, stiff-inplane proprotor 
(designed and constructed by the Bell Helicopter Company and tested in the Ames 
40 by 80 Foot Wind Tunnel in July 1970). The second is a 7.9 m (26 ft) diame- 
ter hingeless, soft-inplane proprotor (designed and constructed by the Boeing 
Vertol Company and tested in the Ames 40 by 80 Foot. Wind Tunnel in August 1972). 
As far as their dynamic characteristics are concerned, the two rotors differ 
primarily in the placement of the nonrotating natural frequencies of the blade 
flap and lag motions. The Bell rotor has a gimballed hub and stiff-inplane 
cantilever blade attachment to the hub, hence vg = 1/rev (nearly, for it does 
have a weak hub spring) and v^ > 1/rev; it also incorporates positive pitch/ 
flap coupling (6^ < 0) to increase the blade flap/lag stability. The Boeing 
rotor has a cantilever or hingeless hub with soft-inplane blade attachment, 
hence vg > 1/rev and v c < 1/vev. The different placement of the blade fre- 
quencies at the opposing extremes of the possible choices, results in quite 
different dynamic characteristics for the two aircraft. 

The rotors are described in references 2 to 5. Table 1 gives the major 
parameters of the rotors and cantilever wings used in the full scale tests (a 
more complete description of the parameters required by the theory is given 
in ref. 1) . The wing frequencies in the theory were matched to the experi- 
mentally measured values by adjusting the spring constants. The typical blade 
frequencies are shown in figures 2 and 3 for the Bell and Boeing rotors, 
respectively. The variation of the Bell lag frequency (fig. 2(b)) with V/fiR 
is due to the collective pitch change. The Boeing rotor blade frequencies 
vary little with collective pitch (V/fiR) since the blade has nearly isotropic 
stiffness at the root. 

The theory was demonstrated to be an accurate representation of proprotor 
and cantilever wing dynamics in reference 1, by comparison with full scale 
wind tunnel test results for the gimballed and hingeless rotors; some of these 
results are also given in reference 6. References 1 and 6 discuss the basic 
(uncontrolled) dynamics of the proprotor and cantilever wing configuration. 
References 1 and 6 also discuss the influence of various elements of the 
theory, including the modelling used for the wing and rotor blade aerodynam- 
ics, the influence of the rotor lag degree of freedom, and the role of the 
rotor rotational speed degree of freedom (windmilling rotor). 


STATE EQUATIONS OF MOTION 

The equations of motion describing the proprotor and wing are of the 

form: 

A 2 Xj ♦ A x Xj + AjjXj* B 0 v ♦ C 0 u 

where 

Xj ■ proprotor degrees of freedom (System State) 
v ■ control components (System control) 
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u = gust components (external disturbance) 

and A 2 , Aj, A 0 , B 0 , C 0 are constant coefficient matrices. 

One or more degrees of freedom may be absent from A 0 x, i.e., a column of A 0 
may be zero (in this case it occurs for the C 0 degree of freedom) . This 
means that there is no spring restraint on these degrees of freedom, making 
the equations for them actually first order, with the results that some of the 
eigenvalues of the second order equation will be zero. It is convenient to 
reduce the system to one of first order. When the substitution 


is made, then 

* 

X2 B *Aj ^^1^2 "^2 ^AqX^ ^ ^2 ^ BqV ^ A2 ^CqU 
and the system becomes 



or, 

x = Ax ♦ B v + Cu (1) 

For the proprotor and cantilever wing system, the state variable vector x is 
comprised of the displacement and velocity o£ each of the nine degrees of 
freedom described above. The input vector v is comprised of the rotor 
collective and cyclic pitch controls (0 O> 6 lc , 6 ls ). Later a wing flap input 
(£p) will also be included in v. For the external disturbance u, only the 
vertical gust component will be considered here. 



SYSTEM PERFORMANCE CRITERIA 


The performance of the system is partly determined from computations of 
its eigenvalues. The system is stable if the real parts of the eigenvalues 
are negative; as the real parts become more negative, the stability of the 
system increases. The damping ratio t (fraction of critical damping) is 


C 



or, on a root locus plot, the cosine of the angle between the negative real 
axis and the line drawn from the eigenvalue to the origin. The system 
becomes more damped as C increases, and is critically damped for C > 1. 

In addition to the stability and damping of the system, the performance 
is evaluated by observing the response of the system to excitation by the 
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vertical component of a gust. The gust used is described by the Von Kdrmdn 
gust spectrum as 

% °w 2l i + 8/3 [1.339 (uL/V) ] 2 
S w (wj = -z-r. k - J ~ 

2 V {1 + [1.339 ( wL / V )] 2 ) n /6 

where u is the frequency (nondimensionalized by £1) ; L the gust correlation 
length (nondimensionalized by R) ; a w the rms value of the gust velocity 
(nondimensionalized by fiR) ; and V is the aircraft speed (nondimensional- 
ized by QR) . L is approximately 1524 n (5000 ft) for high altitude and from 
120 to 150 m (400 to 500 ft) for low altitude. A lower value of L gives a 
stronger high frequency content in the gust spectrum; i.e., the break fre- 
quency, ~ V/L, is larger. For this analysis the L used was 122 m ^400 ft). 
More details of this gust model are given in references 7-10. 

The frequency response of the system to vertical gust is given by 

iut uf \ iut 
xe = H(io)c» G e 


with the transfer function for frequency to. 


H(u) * [-A(A 2 + u 2 I)" 1 C] ♦ i[-u)(A 2 ♦ u 2 I) _1 C] 

where A and C are matrices of (1) . The static response is obtained by 
calculating H(u) for u « 0 [H(0) * -A _1 C]. 

The power spectrum of the response of the state variable B is 
Sg ■ |H(w) | 2 S W , hence the rms response eg is given by 

Og 2 ■ f Sg du ■ T |H(w) | 2 S w du 


(see refs. 7-9). Since Sg is an even function. 



du * 2 J Sgdu 


Computational values for the intestation limits o and • were numerically 
determined to be e" 6 (V/L) and e®(V/L) , respectively, as using more extreme 
limits did not contribute significantly to the value of the integral. 

The ns levels of velocity Og and acceleration og are similarly 
given by 




du 
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du 


Also, the number of peaks, or zero crossings with positive slope, per rotor 
revolution is calculated by 


This equation for N comes from the expression for the expected number of 
crossings of the threshold a per revolution, assuming a Gaussian probability 
distribution, 

6 

with a ■ 0 for threshold level of zero. The ratio of the number of cross- 
ings of the level a to the number of zero crossings gives the fraction of 
cycles in which the response 8 exceeds threshold level o: 

f . e -l/2(a/o 6 )2 

So, to determine the fraction of cycles in which the response 8 exceeds the 
rms level Og, 


- 1/2 

f ■ e ' ■ 60.7% of the cycles. 

Similarly, to find the percent cycles in which the response 8 exceeds n 
multiples of the rms response 


f . ,-" 2 ' 2 

is computed. For example, for n ■ 2, 8 exceeds 2og 13.5% of the cycles; 
for n ■ 3, 8 exceeds 3og 1.1 percent of the cycles. 


OPTIMAL CONTROL PROBLEM 

« 

Given the system (1), the control v is input to the system whose output 
is then measured by the response to an external disturbance u; a system of 
this form is open loop. In a closed loop system, a control $*, based on the 
response of the system is input, or fed back to reduce the response and 
increase the system stability; i.e., the feedback control v* links the out- 
put and the input of the system, closing the loop. The system then becomes 

x ■ Ax ♦ B(v ♦ v*) *(51 (2) 
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The control v* to be determined her. will be an optimal control law, 
optimal in the sense that the maximum displacements of state and control are 
bounded, thus minimizing the measure of performance, or performance index, 


J = J* (x T Qx + v* T Rv*)dt (3) 

o 



The weighting matrices Q and R are somewhat arbitrary, as the magnitudes of 
the elements are n ot important; however, the size of R relative to Q, or 
the 'gain' * /Q/R, determines the level of feedback control. If^ R is large 
relative to Q, i.e., small gain, then to minimize J, control v must be 
small. This allows large values of x and is therefore loose control. If 
R is small relative to Q, i.e., large gain, then control v must be large 
for minimum J, allowing only small values of x and, hence, is tight control. 
Also important are the magnitudes of the elements of Q (and R) with respect 
to each other. If one state degree of freedom is to be controlled more than 
another, then its corresponding element of Q will be larger. Similarly, if 
one control is to be applied more than another, its corresponding element of 
R will be smaller. For this study, Q and R were chosen for simplicity to 
be diagonal matrices. In addition, the elements of Q corresponding to the 
derivative degrees of freedom are set to zero, since attempting to control the 
velocities while already controlling the displacements is redundant; it also 
amounts to a constraint on the frequency content of the state response as 
well as its magnitude, thus introducing the complexity of an additional 
parameter determining Q. Using the maxim>im displacement interpretation, Q 
and R may be represented by 


Q - 


In the presentation of the optimal control results to follow, numbers will 
be given for Q and a gain. For example, the control for a nine degree of 
freedom system may be 

Q » S*.l, 4*1. (equivalent to .1, .1, .1, .1, .1, 1., 1., 1., 1.) 
and GAIN - 10. 

The numbers given for Q are the diagonal elements of the matrix Q corre- 
sponding to the displacement states; all other elements of the Q matrix are 
zero, including the diagonal elements corresponding to the velocities of the 
degrees of freedom. 
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For the matrix R, each control variable is given the weight (GAIN)' 2 ; 
hence R is a diagonal matrix 


(GAIN)' 


where I is the identity matrix. 


Given Q and R, the optimal control problem is to find the control v* 
for the system (2) which minimizes the performance index J, equation (3). 
Details of the formulation of the optimal control problem and its solution are 
given in references 11 and 12. This optimal control law is found to be linear 
feedback of the state: 

v* » -Kx with K » R” 1 B T P 

The patrix P is the solution of the matrix Riccati equation 

PA ♦ A T P - PBR' 1 B T P ♦ Q - 0 

Following references 11 and 12, to solve the Riccati equation for p, the 
matrix 


M . f A - BR ‘ 
-Q -A 


is constructed from the known matrices A, B, Q, and R. M has the property 
that its eigenvalues occur in pairs of equal magnitude but opposite sign. 
With the eigenvectors of M as columns, the matrix 


H ’•] 


is constructed where the columns of T. are the eigenvectors corresponding 
to the eigenvalues of M with negative real parts, and the columns T+ are 
the eigenvectors corresponding to the eigenvalues with positive real parts. 
Since there are equal numbers of the eigenvalues having negative and positive 
real parts, T. and T+ have the same dimensions. T. can be further partioned 


"x_“ 


with X. and Y. having the same number of rows. The solution P is then 


P - Y.X." 1 

Now having sol. ad for v* ■ -Kx, the do* d loop system (2) becomes 
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Ax ♦ Bv - BKx ♦ Cu 


or 


* 

x = (A - BKJ x ♦ Bv ♦ Cu . 

The matrix (A - BK) then replaces the A of (1) and the performance of the 
closed loop system may be evaluated and compared to the open loop performance. 


RESULTS AND DISCUSSION 

The optimal control theory methods will be used to design linear state 
variable feedback for two proprotor/canti lever wing designs -- a gimballed, 
stiff-inplane rotor and a hingeless, soft-inplane rotor. The rotors are 
assumed to be in windmilling operation. That should be a good model for the 
dynamics considered here since the enginc/drive-train/govemor dynamics will 
nave a long time constant compared to the proprotor and wing modes. For the 
frequency range of interest here the rotor should behave as if windmilling 
even with a rotor-speed governor operating. This assumption may be less valid 
for the system longitudinal dynamics, where the rotor rotational speed degree 
of freedom has a particularly important role. Therefore this analysis is 
concerned primarily with the lateral/vertical dynamic of the proprotor and 
cantilever wing system — for example, the response of the wing vertical 
bending (q.), tip path plane tilt iB Jc and B ls ), and cyclic lag (c lc and ; ls ) 
to vertical gusts. It should be noted however that the optimal control 
solution includes reedback and control of the longitudinal dynamics as well; 
the vertical gust input enters only in examining the system performance, not 
at all in the optimal control solution. 

The primary criteria for the system performance will be the root mean 
square (rms) response of the rotor motion (og and o^) , the rms acceleration of 
the wing motion (o" , n- , and o-) , and the damping ratio of the wing modes 

Hj ^2 * 

(q,» Q 2 * P)* ^ irst 15 8 measure of the rotor motion for an articulated 

rotor, or of the 1/rev blade loads for a hingeless rotor; the second is a 
measure of the airframe vibration due to the gust; and the last is a measure 
of the overall system dynamic stability. Normally the rms response and 
acceleration will be given as a ratio of closed loop to open loop values; 
occasionally the absolute magnitudes of the motion will also be examined. For 
the rotor motion, the rms value of the amplitude of the cyclic flap and lag 
motion is given — the square root of the sum of the squares of the 1C and IS 
degrees of freedom (|e| ■ (6 1C 2 ♦ S ls 2 ) 1 ^ 2 end |c| • (; 1C 2 ♦ Ci« 2 W 2 - In 
addition to the damping ratio of the wing modes, complete root loci diagrams 
for the system will be given. 

First, the optimal control results will be examined for a standard case 
to determine the improvement in the dynamics possible. Then the results will 
be examined for cases when the controller is designed considering only the 
wing or only the rotor. Next the effect of adding a wing flap to the avail- 
able controls will be considered. Finally, the effect of cruise speed on the 
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control system design will be investigated, considering the performance with 
an optimal controller designed for each speed, and the performance with a 
controller designed for a single speed but operated for all speeds. 

Open Loop Transfer Functions 

Figures 4 and 5 show the open loop transfer functions for the Bell and 
Boeing rotors respectively (at V/fiR = 0.7, windmilling rotor). The magnitude 
of the response of each of the nine degrees of freedom to a vertical gust 
angle is given. Also shown are the transfer functions multiplied by 

/S w /S w (0) (where S w is the gust spectrum, given above) . The square of 

|H| • is the spectrum of the response of the system to the gust, and its 

integral is the rms response. Note the prominence of the wing vertical bending 
mode resonance in the transfer function, even when multiplied by the gust 
spectrum. This implies that the wing and rotor dynamics have a significant 
role in the res ponse, and hence in the control required. Figure 6 shows the 
gust spectrum, /S w /S w (0) for V/J2R = 0.7. The gust spectrrm has a very low 
comer frequency (around 0.025/rev) even at this high speed (i.e., the gust 
wavelength is very long) . However, the square root of the gust spectrum only 
falls off as S* 1 / 2 ~ w 5 / 6 . Hence the high frequency response — the reso- 
nances in the transfer function — contributes significantly to the gust 
response. 

Gimballed, Stiff- Inplane Rotor 

For the gimballed stiff-inplane rotor (Bell) and wing, Q = 5*.l, 4*1 and 
V/OR =0.7 was taken as the standard case. The values of Q were chosen so 
that the rms acceleration due to vertical gust achieved a fractional reduction 
that was roughly the same for both the wing and rotor degrees of freedom. The 
gain was varied from 0 to 40 with most of the effect occurring by 10 and 
leveling off from there. Figure 7 shows the results of this case. Part 7(a) 
shows that a substantial reduction in Tins flap and lag motion is possible by 
applying suitable control. The reduction in the rms flap response is from 
0.74 deg/deg at zero gain to around 0.17 deg/ deg at high gain, indicating a 
77 percent reduction in flap amplitude; so then, for a 6 m/sec (20 fps) rms 
gust (strong turbulence), the reduction in rms flap response is from 2 deg to 
0.46 deg. Figure 7(b) shows that the number of peaks/revolution for B and c 
remain about the same order as the gain increases; however, there is some 
increase due to o x being reduced faster than o x . 

For wing vertical bending qj, the dominant wing mode, figure 7(c) shows 
that the rms response achieves about a 45 percent reduction which for a 6 m/ sec 
(20 fps) gust would correspond to reducing the rms vertical motion at the wing 
tip from about 4 cm (1.5 in.) to about 2 cm (0.8 in.). Also shown is the num- 
ber of peaks/revolution (also in cycles/hour) and its corresponding reduction. 

A result of these response and peaks/rev reductions is a decrease in the cumu- 
lative wing fatigue. In figure 7(d) is shown a 55 percent reduction in the 
rms acceleration of q 2 which for a 6 m/sec (20 fps) gust reduces the rms 
vertical acceleration at the wing tip from about 1.25 g to about 0.57 g. 
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Figure 7(e) shows that both the wing and rotor rms acceleration due to 
vertical gust reduces significantly (55 percent to 75 percent) as gain 
increases. However, while the rotor motion (°g» <*») continuously decreases 
with increase in gain, the reduction of the wing motion tends to reach an 
asymptotic limit. This implies that for a given Q, there is a limit to the 
reduction in wing motion possible, possibly due to the fact that the only 
controls are to the rotor and the ratio of wing motion to rotor motion is 
fixed (because of fixed Q) . Figure 7 (f) shows that the rms response exhibits 
similar behavior to the acceleration with a 45 percent to 75 percent reduction 
achieved. In figure 7(g) there are again significant reductions in the steady 
state response as gain increases; however, in many cases the reduction is much 
less than for the rms response or acceleration. In all cases, the static 
response is not a good measure of rms response to gust, implying that the 
higher frequency dynamics have a large role in the gust response, especially 
the wing vertical bending resonance peak. 

The damping ratios of the wing modes shown in figure 7(h) indicate a 
substantial increase; e.g., the wing vertical bending mode (qj) damping 
increases from about 3.5 percent to about 12 percent of critical damping. 
Although whe wing modes in particular show a leveling off at higher gain, the 
limit on wing control may not be a critical factor, ai the damping increase 
achieved is substantial and probably all that would be required. However, 
this leveling off at high gain may be associated with the type of control and 
the criterion used, i.e., state variable control and the quadratic performance 
index. 


The feedback gain matrix K in the expression 
with gain = 10 is shown in table 2. 


v = -Kx for this case 


It was hoped that the matrix would show that measurements of only a few 
quantities and feedback to only a few inputs would dominate the control, i.e., 
many of the elements of K would be insignificant compared to a few dominant 
ones; for then the matrix would suggest a simpler control scheme (smaller 
matrix). Instead it is found that many of the 51 elements of K are of the 
same order, and few control elements could be considered unnecessary. Part 
of this complexity is due to the fact that the whole system is being con- 
trolled, not just the dynamics due to the response to vertical gust, which is 
not a design criterion but a measure of performance. The gain matrix is not 
without some pattern however; namely, there is some characteristic decoupling 
into lateral/vertical and longitudinal systems: 0 O , collective control, is 

obtained mainly from measurements of 6 0 , q 2 , £ 0 , g 0 , q 2 , the longitudinal 
degrees of freedom; and 6 lc /9 1 e, cyclic control is obtained mainly from 
feedback of i n is , q, p, 6, As. c /£,_, qj, p, the vertical/lateral 
degrees of freedom. Even with this decoupling, the feedback is still quite 
involved however, but it does suggest the possibility of designing a con- 
troller for a reduced system; e.g., analyzing and designing a control system 
for the longitudinal dynamics on the basis of just the q 2 , 6 , and £ 0 degrees 
of freedom. A controller of this type must always be checked in operation 
with the complete system dynamics however. 
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The root locus in figure 8 further illustrates the effects discussed. 

There is an increase in the damping of the wing modes, leveling off at high 
gain. Also there is a continual large increase in the damping of the rotor 
lag modes but the flap modes show little significant change, though the flap 
and lag motion are highly coupled by the rotor aerodynamics due to the high 
inflow ratio. Because of this coupling, the labeling of the 8-1 and c-1 
modes is somewhat ambiguous (on this plot as well as the other root loci) with 
the further coupling introduced at high gain. 

Figure 9 shows the response (transfer functions) of the degrees of free- 
dom to a vertical gust for gains of 0, 2, 10, and 40. There is a general 
reduction in the magnitude of the response as gain increases; in particular, 
there is a reduction in the wing vertical bending resonance peak in each 
degree of freedom. Note also the reduction of the response at high frequencies 
(around 2/rev) with very high gain. 

From the case of Q = 5*.l, 4*1., it is concluded that simultaneous, 
significant reductions in rotor and wing response with an accompanying 
increase in dynamic stability of the system is achievable. The feedback con- 
trol required is not simple as many elements of the state variable feedback 
gain matrix are of the same order. With only the rotor controls as inputs to 
the system, there appears to be a limit to both the reduction in wing motion 
and the increase in wing damping possible (due possibly to the control system 
design technique used, not the proprotor configuration) for although the rotor 
response continues to decrease at high gain, the wing motion is asymptotic 
to nonzero levels. However, the level of reduction of the wing motion possible 
using just rotor controls may be satisfactory. 

The next case examined was for V/flR =0.7 and with Q = 6*1, 3*0, 
i.e., only the rotor degrees of freedom have constraints on their motion so 
that the wing motion does not contribute to the quadratic performance index. 

In figure 10(a) the rms acceleration shows very little reduction is achieved 
in wing response (q 2 , q,, p) using this performance ind x, especially for wing 
vertical bending (q.) which has a reduction of only 5 pv cent to 10 percent. 

In contrast, the reauction of the rotor flap and lag acceleration is greater 
than for the previous case with the wing motion constrained also; there is 
about 75 percent reduction at high gain compared to about 65 percent for the 
previous case. The rms response in figure 10(b) has characteristics similar 
to those of the acceleration. Wing motion shows little reduction, e.g., qj 
is reduced only 10 percent to 15 percent. However, q 2 actually increases in 
rms response due to coupling with the 8 0 and c 0 degrees of freedom (but 
since the absolute level of the q* 2 response is still quite small, this is no 
real problem). With a slightly greater weight on £ 0 (smaller corresponding 
value of Q), the effect could be eliminated. The rotor flap and lag rms 
response is reduced significantly by elimination of the constraint on the wing 
motion, with about an 85 percent reduction at high gain compared to about 
75 percent for the previous case. 

The variation of the damping ratio of the wing modes with gain in 
figure 10(c) especially shows the effect of eliminating the constraint on the 
wing motion. The damping of the q 2 and p (wing chordwise bending and wing 
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torsion) modes change very little and the damping of the qj (wing vertical 
bending) mode actually decreases, from about 3.5 percent for the open (gain = 0) 
level down to about 2.5 percent at higher gain, as compared to an increase to 
about 12 percent reduction achieved in the previous case. This decrease in 
the damping of the wing vertical bending mode is very serious as the dynamic 
stability of the wing modes is a critical factor in proprotor aircraft design, 
and implies that a control system designed considering the rotor only will 
not be satisfactory. 

The gain matrix, shown in table 3, as compared to the previous gain matrix 
(table 1), shows a significant decrease in feedback of wing motion; e.g., the 
reduction in the and qj elements. This gain is effectively three times 
the magnitude of the previous gain matrix (because of the greater weighting 
of the rotor degrees of freedom) and some of the differences are obscured. 

In figure 10(d), the root locus shows a variation of the rotor modes 
(8 ± 1, c ± 1, 8) with gain that is almost identical to the previous case. 

But with no constraint, the wing modes (q^, q 2 , p) change very little with gain 
and do not show the significant increase m damping as with the first Q. 

The next case considered used V/ftR = 0.7 and Q = 6*0, 3*1, i.e., con- 
straints are on the wing degrees of freedom only, with no contribution of the 
rotor motion to the performance index. Figure 11(a) shows that there is a 
greater reduction in the wing rms response in this case compared to the first 
case, with about a 65 percent reduction in qj compared to about 45 percent 
for the first Q. Here especially is seen the effect of relaxing the con- 
straint on the rotor motion; there is a 40 percent reduction of the rotor 
motion response, much less than the 75 percent reduction for the first Q; 
also at higher gain, the rotor response then increases again and even becomes 
significantly greater than the open loop (gain = 0) level for gain above 20. 
Figure 11(b) indicates that a somewhat greater reduction in wing vibration is 
achieved by relaxing the constraint on the rotor motion, a 65 percent reduc- 
tion in qj compared to 55 percent for the first Q. However, the cost of 
this increased control of wing motion is a much smaller reduction in the rotor 
motion. 

Figure 11(c) shows that there is a tremendous increase of the damping of 
the wing modes with increasing gain with no weight placed on the rotor 
motion, although initially the increase is about the same as with the first 
Q. The qj damping increased from 3.5 percent of critical damping to around 
35 percent at high gain, compared to j 12 percent increase achieved with the 
first Q. 

The gain matrix in table 4 shows that there is almost no feedback of the 
rotor motion at all, just for the wing degrees of freedom, q 1( q 2 , and p. 

The feedback law is then primarily 

e is - -Mi - M 

e o " M 2 
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(with some displacement of wing motion feedback as well) which then explains 
the large wing mode damping increase. This is as expected for increasing the 
wing damping, but with no constraint on the rotor motion, the rotor response 
increases with the wing damping. 

Figure 11(d) shows that the wing modes (qj, q„, p) show dramatic, large 
increases in damping as gain increases, as compared with the first Q. With 
no constraint on the rotor modes, their roots do not change much, with the 
flap mode roots (0 ± 1) actually showing some decrease in damping which, 
however, slows down at high gain. 

It is seen that relaxing the constraint on the wing motion achieves a 
larger reduction in the rotor response at the cost of less control of the wing 
motion (the wing vertical bending mode damping actually decreased) and vice 
versa for relaxation of. the constraint on rotor motion (flap rms response is 
actually increased) . The gains achieved by relaxing the constraint on part 
of the system are in fact much less than the cost; modest improvement in the 
response of part of the system is possible but only with a resulting deteriora- 
tion in the dynamic characteristics of the rest of the system. These results 
may be taken as limiting cases of possible apportioning of the performance 
index weight between rotor and wing. The optimal control procedure here pro- 
vides the means to design a control system considering the entire system, with 
the first Q selected being a good intermediate case between the two 
extremes. With a proper apportioning of weight in the quadratic performance 
index to the wing and rotor, it is possible to simultaneously improve the 
dynamic behavior of both the rotor and wing, achieving nearly the levels 
obtained when designing a controller for the rotor or wing alone. So, to 
improve all the dynamic characteristics of the proprotor aircraft, it is 
necessary to consider the complete system, both wing and rotor, although the 
feedback control that results from this analysis is not simple. 

T 'ie next case considered used Q = 5*.l, 4*1., V/ftR = 0.7, and a gain 
range of 0 to 40, as in the first case Furthermore, added to the controls 
was a wing trailing edge flap. A 30 percent chord flap was used, extending 
over the outer 50 percent span of the wing. The effect on the aerodynamic 
lift and pitching moment on the wing were calculated by strip theory analysis. 
The control vector for this case now has 4 elements: 0 Q , 6 lc , 6 1$ , 6p. 

Figure 12(a) shows that a substantially greater reduction in wing verti- 
cal bending rms acceleration is achieved by the addition of the wing flap 
control, about 85 percent reduction as compared to about 55 percent for the 
fir- . case (no flap control) , with. less of a tendency for the wing response to 
livel off at high gain. Comparing the reductions in the case of a 6 m/sec 
(20 fps) rms gust, the rms vertical acceleration at the tip decreases from 
0.57 g in the first case to about 0.18 g with wing flap control (cf. 1.25 g 
uncontrolled). This reduction in wing response is achieved with even less 
rotor response, reduced by about 70 percent here as compared to 65 percent 
reduction for the first case. 

The rms response in figure 12(b) again shows .. reduction in wing and 
rotor response of about 85 percent at high gain as compared to about 75 percent 
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for the first case (only rotor control). For moderate gain levels, about 
95 percent reduction is achieved in flap response; applied to a 6 m/sec (20 fps) 
gust, that gives Og = 0.1 degrees as compared to og = 2 degrees *n the 
uncontrolled case. 

In figure 12(c) the wing mode damping shows dramatically the effect of 
including the wing flap control. At high gain the wing vertical bending mode 
dampi 'g is about 50% of critical damping as compared to about 12 percent 
without flap control, and 3.6 percent uncontrolled. The torsion (p) damping 
shows a similar increase; however, the wing chordwise bending mode is not 
affected, as would be expected since the flap gives primarily lift and torsion 
movement on the wing. 

The feedback to the controls shown in table 5 is similar to the case 
without the flap control, except that there is less feedback of the wing 
motion, as that feedback now goes to the flap. The flap control is primarily 

6 F = -K 1 q 1 

which directly increases the wing vertical bending damping. So, the flap 
takes over the primary control of the wing motion (except the chordwise bend- 
ing which is still controlled by the rotor collective) with little feedback of 
the rotor motion. 

In figure 12(d) the root locus is nearly identical to the case without 
flap control, except that the wing vertical bending and torsion modes have 
greatly increased damping. 

It has been shown that the wing flap is very powerful for the control of 
proprotor/ cantilever wing motion. The forces due to the flap deflection act 
directly on the wing to control wing motion, especially the important wing 
vertical bending motion. Simple feedback of the modal velocities allows a 
direct increase of damping to very high levels (on the order of 50% critical 
in this example). By taking over a major role in controlling the wing, use of 
the flap then allows more control of the rotor for a given level of rotor 
control deflection (i.e., gain), thus allowing a greater reduction in rotor 
motion as well as wing motion. Therefore, it is probably desirable to make 
use of wing aerodynamic control surfaces in the feedback control of proprotor 
aircraft dynamics. However, the effectiveness of such control surfaces 
depends on the viscous flow, compressibility, unsteady aerodynamics, three- 
dimensional aerodynamics; structural and inertial limitations on flap oscilla- 
tion frequency; and other effects.' Hence, such controllers will certainly 
require experimental verification. 

The last case examined for the Bell proprotor used Q = 5*.l, 4*1., and 
gain = 10 while varying the velocity from 100 to 400 knots. There are two 
reasons for considering this case; the first is to examine the variation in 
response of the rotor and wing due to changing cruise speed, both open loop 
and optimally controlled; secondly, to determine if it is necessary to vary 
the controller design with speed in order to always achieve the desired 
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response, i.e., to examine the performance of the controller operating at a 
speed other than the one at which it was designed, hence the expression, 

"off design." For this latter case, an optimal controller is designed for a 
particular speed and the response using that controller at several other 
speeds is calculated. Compared to that response is the response calculated 
using the optimal controller for each speed. For this analysis the off 
design performance of a 250 knot controller will be examined. 

Figure 13(a) shows the rms acceleration using the optimum controller for 
each speed. There is little variation in the percent reduction of the 
response for a given gain with change in speed, with the exception of the wing 
vertical bending response which reaches a peak at 200 knots; at that speed 
there is a resonance of the wing and rotor motion, specifically, between the 
wing vertical bending mode (qj) and the rotor low frequency lag (5-1) mode. 
Away from that resonance, which produces substantial coupling of the wing and 
rotor motion, a greater percent reduction of the q : motion is possible for 

a given gain. 

Figure 13(b) shows the rms acceleration due to using the 250 knot con- 
troller at each speed instead of the optimal controller (fig. 13(a)); the 
response is nearly the same for both cases. 

The rms response using the optimal controller at each speed in fig- 
ure 13(c) exhibits similar behavior to the acceleration. Figure 13(d) shows 
that using the 250 knot controller similarly produces nearly the same response 
as in figure 13(c), except for the chordwise bending (q 2 ) response (though 
the actual magnitude of the response of that mode to vertical gust is small). 

Figure 13(e) shows improvement in the wing modes damping levels for all 
speeds using the optimal controller at each speed. Using the 250 knot con- 
troller as shown in figure 13(f), gives about the same results. In both 
cases (optimal controller and 250 knot controller), the controller is capable 
of substantially increasing the wing vertical bending mode damping (and other 
modes) even at high speed, while in the uncontrolled system, the damping 
decreases as speed increases due to rotor high inflow aerodynamics. 

The gain matrices for the optimal controller at 100, 250, and 400 knots 
in tables 6, 7, and 8 respectively, show very little variation in either 
absolute or relative magnitudes of the elements, accounting for the success 
of the 250 knot controller over the entire speed range. 

Figure 13(g) shows the root locus plot for both the uncontrolled and the 
optimally controlled system for varying speed and indicates general improve- 
ment of the system stability at all speeds with the optimal controller. In 
figure 13(h) is shown the root locus for the uncontrolled system and the sys- 
tem using the 250 knot controller for varying speed, demonstrating nearly the 
same results as with the optimal controller at each speed for both wing and 
rotor modes. 

Figure 13 (i) shows the actual magnitudes of the response of the flap and 
lag motion in deg/fps of gust while varying speed for a) the uncontrolled 
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system, b) using the optimal controller at each speed, and c) using the 
250 knot controller at each speed. The uncontrolled system shows a tremendous 
increase in rotor response with increasing speed due to the high inflow aero- 
dynamics. While the controlled system also shows an increase in actual magni- 
tude of the response with increasing speed, it does so at a much slower rate. 
Applied to a 6 m/sec (20 fps) gust, the open loop flap and lag rms response at 
400 knots is about 2.5 degrees and is reduced to about 0.5 degrees using a 
controller (at gain - 10). 

Figure 13 (j) is similar to 13 (i) except that it shows the magnitude of 
the vertical acceleration of the wing tip in g/fps of gust. The uncontrolled 
system shows a great increase in qj with speed, even faster than the rotor 
motion increases, while the increase in the controlled system response is much 
slower (recall that the fractional closed loop/open loop level decreases as 
speed increases from 200 to 400 knots, figs. 13(a) and (b)). For a 6 m/sec 
(20 fps) gust then, the rms vertical acceleration at the wing tip for 400 knots 
is about 2.9 g for the uncontrolled system and is reduced to about 0.9 g with 
a controlled system (gain = 10). 

It can be concluded that a substantial reduction in system response and 
improvement in dynamic stability over the entire cruise speed range of the 
aircraft may be achieved with a properly designed controller. The results and 
conclusions in the rest of this study, for the V/f<R = 0.7 case may then be 
applied to the entire speed range. It has been shown that a controller 
designed for a particular cruise speed (250 knots here) has excellent perfor- 
mance characteristics at all speeds, nearly the same, in fact, at off-design 
speeds as the performance with the optimal controller for that speed. 

Hingeless, Soft- Inplane Rotor 

The behavior of the controlled system and the conclusions based on this 
optimal control analysis for this rotor are much the same as for the gim- 
balled rotor. This discussion therefore will concentrate on those points 
particular to the case of the Boeing rotor. The differences between the two 
rotors are primarily the flap and lag frequencies, the Boeing rotor having 
cantilever, soft-inplane blades, while the Bell rotor has a gimballed hub 
with stiff-inplane blades. The hub moment capability of the Boeing rotor does 
have a significant effect on the dynamics, both open loop and closed loop. 

For the hingeless, soft- inplane rotor and wing, Q = 5*. 2 , 10, 1 , 10, 1 
and V/S1R = 0.7 was taken as the standard case. The gain was varied frv..n 
0 to 40 with most of the effect occurring by 10 and leveling off from there. 
Figure 14 shows the results of this case. Figure 14(a) shows that the flap and 
lag rms response achieve a substantial reduction over the closed loop values. 
The reduction in the rms flap response is from 0.44 deg/deg at zero gain to 
around 0.15 deg/deg at high gain, about a 65 percent reduction in flap ampli- 
tude. Then for a 6 m/sec (20 fps) rms gust, og (rms flap response) is reduced 
from 1.4 degrees to 0.48 degrees and (rms lag response) is reduced from 
3.1 degrees for the uncontrolled system to about 0.28 degrees. Figure 14(b) 
shows that the number of peaks/revolution remain about the same order as gain 
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increases, but there is some increase due to o x decreasing faster than o*. 

Since the number of peaks/rev in the nonrotating frame is low (fig. 14(b)), it 
follows that the peaks/rev in the rotating frame is dominated by the blade 
once-per-rev motion, for both the uncontrolled and controlled systems. 

Figures 14(c) and (d) show a significant reduction in rms response and 
acceleration of the wing vertical bending motion as gain increases, with about 
60 percent reduction achieved in the rms response and about 70 percent in the 
rms acceleration. For a 6 m/sec (20 fps) gust, this corresponds to reducing 
the rms vertical displacement at the wing tip from 5 cm (2.1 in.) to 2 cm 
(0.8 in.), and the acceleration at the tip from about 1.06 g to about 0.35 g. 

Also shown is the number of peaks/revolution (and cycles/hour) and the corre- 
sponding reduction. 

The ms acceleration shown in figure 14(e) achieves a significant reduc- 
tion in both rotor and wing response to vertical gust. Again, it is seen 
that while the rotor motion continues to decrease with gain, the wing motion 
reduction tends to level off at nonzero levels for higher gain. A reduction 
of 65 percent to 80 percent is achieved (ignoring the q 2 response which is 
always small for vertical gust). Figure 14(f) similarly indicates a 65 percent 
to 90 percent reduction achieved in the rms response. The steady state 
response in figure 14(g) is reduced as gain is increased, but is not at all 
similar, indicating the important role of the high frequency dynamics in the 
gust response. 

Figure 14(h) shows that there is a substantial increase in the wing 
modes damping, e.g., wing vertical bending (q x ) damping is increased from 
about 1.2 percent of critical damping to about 6.7 percent. The open loop 
(gain *= 0, or uncontrolled) low damping level for the q t mode is due to the 
air resonance behavior of this soft-inplane rotor/cantilever wing configura- 
tion, but this mode can be significantly stabilized while simultaneously 
reducing the rotor response by use of the appropriate feedback control. 

The feedback gain matrix in table 9 shows the increased importance of the 
rotor lag mode velocity and displacement in the lateral/vertical dynamics 
(the feedback to the controls 8 lc /©, s ) due to air resonance behavior; i.e., 
it is the coupling of q 2 with c lc /5 ls which produces low open loop 
(uncontrolled) stability, so the controller feeds back lag motion as well as 
wing vertical bending to stabilize this effect. The rms response and acceler- 
ation reflect this control, as there is a very large reduction in rotor lag 
motion. 

The root locus shown in figure 15 indicates an increase in the system 
stability evident in the general movement of the roots to the left. The 6-1 
and (-1 modes are highly coupled and, since the q } mode also has about the 
same frequency, it couples with the other two (more so in other root loci to i 

follow). Consequently, the labeling of these modes at high gain is somewhat 
arbitrary. The frequency response to vertical gust in figure 16 shows a 
general reduction in the magnitude of the response as gain increases with a 
large decrease in the wing vertical bending mode peak and decrease? also in 
the high frequency peaks. 
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These results indicate that simultaneous, significant reductions in 
rotor and wing response and an increase in the stability of the system are 
possible. 

The next case considered used V/flR = 0.7 with Q = 6*1., 3*0., placing 
a motion constraint on the rotor degrees of freedom only. The rms accelera- 
tion in figure 17(a) shows that little reduction in the wing motion is 
achieved; particularly, the wing vertical bending is reduced only about 
30 percent. There is, however, a greater reduction in the rotor motion in 
this case than when the constraints were applied to the wing degrees of free- 
dom also, about 75 percent reduction in this case as compared to 55 percent 
in the first case. Figure 17(b) shows only about 30 percent reduction in the 
response of qj, but about 75 percent reduction in flap response compared to 
65 percent in the first case. The wing damping in figure 17(c) shows little 
change with increasing gain. q 1 damping remains at the low level of 1.5 per- 
cent critical. The gain matrix shown in table 10 shows a reduction in feed- 
back of the wing motion, especially ^ and q Jf compared to the first Q. 

Figure 17(d) shows that the rotor roots vary nearly the same as for the first 
case while the wing roots change very little. 

Next was considered the case of Q * 6*0., 3*1., having constraints on the 
wing motion only, and V/flR * 0.7. The rms acceleration in figure 18(a) shows 
that a somewhat greater reduction in wing vibration is achieved, about 
75 percent reduction in compared to 65 percent for the first case. But 
the cost of this reduction is that the rotor response is not reduced as much 
as when there are constraints on the rotor motion also. Figure 18(b) illus- 
trates especially the cost of increased control of the wing motion. The rotor 
flap response decreased by about 55 percent, compared to 65 percent reduction 
in the case of the first Q; the rotor lag response reaches its maximum 
reduction of about 40 percent at a gain of 5 as compared to a previous 90 per- 
cent reduction, and then increases again at higher gain, becoming greater than 
the open loop value for gain greater than 25. The wing damping in figure 18(c) 
experiences a tremendous increase with no constraint on the rotor motion. The 
q, damping increases from around 1.2 percent of critical damping open loop to 
about 45 percent at high gain, compared to about 7 percent achieved with the 
first Q. The gain matrix in table 11 shows almost no feedback of the rotor 
motion. The feedback law is primarily 

e ls . -k^ - k 2 * 

»0 * Ms 

with some feedback of wing displacement as well. In figure 18(d), the root 
locus indicates a large increase in the damping of the wing modes. Because of 
their similar frequencies, there is significant coupling of the 8-1 , c-1, and 
qj modes at higher gain. 

It has been shown that relaxation of the constraint on the wing notion 
achieves a larger reduction in the response of the rotor modes, but with a 
resulting decrease in the control of the wing motion and vice versa for 
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relaxing the constraint on the rotor motion. Thus it is necessary to consider 
both the rotor and wing motion when designing a controller so as to simultan- 
eously greatly reduce the response to vertical gust and increase the stability 
of the system. 

In the next case a wing trailing edge flap is added to the controls. As 
for the first case, V/ft = 0.7, Q = 5*. 2, 10, 1, 10, 1 and the gain ranges 
from 0 to 40. Figure 19(a) shows that, with this added flap control, there is 
a substantially greater reduction in wing vertical bending motion, 90 percent 
reduction at high gain compared to about 70 percent without using the flap con- 
trol. Applied to a 6 m/sec (20 fps) rms gust, the rms vertical acceleration 
at the tip is decreased from 0.35 g with no flap control to 0.11 g with flap 
control (cf. 1.05 g uncontrolled). Shown in figure 19(b) is the reduction in 
wing response with a simultaneous reduction in rotor motion by using the wing 
flap control. The reduction is around 85 percent compared to 65 percent with- 
out this control. So, using this control for a 6 m/sec (20 fps) rms gust 
gives eg = 0.2 degrees versus 1.4 degrees for the uncontrolled case. 

Shown in figure 19(c), wing vertical bending and torsion modes damping 
experience a dramatic increase (chordwise bending mode damping unchanged); 
about 55 percent of critical damping is achieved for the wing vertical bending 
mode at high gain compared to only about 7 percent without wing flap control. 
The gain matrix in table 12 shows that little of the rotor motion is fed back 
to the wing flap control, but there is a transfer of the wing q, and p 
feedback from the rotor to the flap control. The root locus in Figure 19(d) 
is nearly identical to the root locus for the case without using the flap, 
except the greatly increased damping of the wing vertical bending and torsion 
modes. 

It is seen that the wing flap is very effective in controllirg the 
proprotor/canti lever wing motion, especially the wing vertical bending motion, 
and allows more control over the rotor for a given gain level, and thus a 
greater reduction in rotor motion, as well as wing motion. 

The last case examined for the Boeing proprotor used Q ■ 5*.l, 4*1., and 
gain » 10 while varying the velocity from 100 to 400 knots. The variation of 
the response with speed was examined, using the optimal controller at each 
speed and using also a 250 knot controller (optimum at 250 knots) at each 
speed. In figure 20(a), there is little influence of speed on the effective- 
ness of the optimal controller. There is some increase in the percent reduc- 
tion of the response at low speed where there exists a resonance of this 
mode with the rotor low frequency -lag mode (air resonance behavior) . The 
response with the 250 knot controller shown in figure 20(b) is much the same 
as with the optimal controller for each speed, except for the wing chordwise 
bending which has only a small response to vertical gust anyway, and a some- 
what increased flapping response at low speed. Figure 20(c) illustrates that 
the rms response behavior is similar to that of the acceleration when using 
the optimal controller at each speed Using the 2S0 knot controller, it is 
shown in figure 20(d) that the response is nearly identical as when using the 
optimal controller for each speed. An exception is the somewhat increased 
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flapping response at low speed as compaied to the optimal controller results, 
but the actual magnitude of the response is low at that speed. 

Shown in figures 20(e) and (f) is a substantial and uniform improvement 
of the wing mode stability over the entire speed range, with nearly the same 
performance using the 250 knot controller for all speeds. In tables 13, 14, 
and 15 are the gain matrices for the optimal controller at 100, 2S0, and 
400 knots respectively. They show little variation, though somewhat more than 
for the Bell rotor, accounting for a somewhat greater variation in performance 
especially at low speed. Figures 20(g) and (h) show the root loci for a) the 
optimal controller at each speed and b) the 250 knot controller. All modes 
show a consistent and substantial improvement, especially the rotor modes with 
the 250 knot controller. 

Figure 20 (i) shows the magnitude of the response of the flap and lag 
motion in deg/fps of gust with variation in speed for a) no controller, b) the 
optimal controller, and c) the 250 knot controller. There is a tremendous 
increase of the response with speed for the uncontrolled system; and while 
there is some increase for the controlled system, it is much slower. Com- 
parison of the optimally controlled system response with the 250 knot con- 
trolled system response shows little difference, indicating that the differ- 
ences at low speed are unimportant. If the results between the uncontrolled 
and controlled systems are compared for a 6 m/sec (20 fps) rms gust, the open 
loop response at 400 knots is about 4.5 degrees for lag and 2.2 degrees for 
flap, and is reduced to about 0.4 degrees for lag and 0.8 for flap closed loop 
(for gain ■ 10). Figure 20(j) is similar to figure 20(i) except that it shows 
the magnitude of the rms acceleration in g/fps of gust at the wing tip, and it 
exhibits similar behavior; the uncontrolled response shows a great increase 
with speed, with a much slower increase using a controller. The 250 knot con- 
troller performance is nearly identical to that with the optimal controller at 
each speed. Applying the results to a 6 m/sec (20 fps) gust, the rms vertical 
acceleration at the wing tip at 400 knots is about 1.7 g uncontrolled and is 
reduced to about 0.6 g with a controlled system at gain ■ 10. 

It can be concluded that a properly designed controller is effective and 
greatly improves the system dynamic performance over the entire cruise speed 
range. Also a controller designed for a particular cruise speed (2S0 knots 
in this case) gives excellent performance, nearly the same as the performance 
obtained using the optimal controller at each speed. 


CONCLUDING REMARKS 

The following conclusions are reached from this investigation, based on 
the proprotor/cantilever wing model and the optimal control theory design 
technique. 

a) Significant and simultaneous reduction in the rotor response, reduc- 
tion in the wing response, and increase in the dynamic stability of the sys- 
tem is achievable. However, the optimal state variable feedback control 
required is not simple. It should be noted that these results are based on 
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feedback of the state variables. In practice the state variables must be 
estimated from some measurement of the system motion, and so measurement noise 
and inaccuracies in the mathematical model of the system will reduce the con- 
trol system performance below the optimum achievable. 

b) In the design of control systems to improve the dynamic character- 
istics of the proprotor aircraft, it is necessary to consider the complete 
system. Considering just the rotor or just the wing dynamics will not produce 
a successful design. 

c) Using a wing flap in addition to the rotor inputs is very powerful 
for control of the proprotor/canti lever wing system. The wing flap acts 
directly to control the wing vertical bending motion, the critical element of 
the dynamics, and so allows very high levels of wing damping to be achieved 
while freeing the rotor inputs for more control of the rotor motion. 

d) A control system designed for a particular cruise speed (250 knots 
here) has excellent performance at all forward speeds, nearly the same in fact 
at off-design speeds as the performance with the optimal controller for that 
velocity. 
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TABLE 1.- DESCRIPTION OF THE FULL-SCALE PROP ROTO PS (AS TESTED IN THE AMES 

40- BY 80-FT WIND TUNNEL) . 



Bell 

Boeing 

Rotor 



Type 

gimballed, 

stiff-inplane 

hingeless , soft- 
inplane 

Number cf blades 

1 

3 

3 

Radius , R 

3.81 m 

1 

3.96 m 

Lock number, y 

3.83 

4. 04 

Solidity ratio 

0* 089 

0.115 

Pitch/flap 
coupling, 63 

-15 deg 

0 

Rotor rotation 
direction , on 
right wing 

clockwise 

i 

j 

counterclockwise 

Tip speed, ftR 
(cruise mode) 

183 m/sec 

160 m/sec 

Rotation speed, ! 

(cruise mode) 

458 rpm 

386 rpm 

' 

Wing 



Semispan, y w / R 

1.333 

1.281 

Mast height, h/R 

0.342 

0.354 

Typical frequencies 
Vertical bending 

3.2 Hz 0.42/rev 

2.3 Hz 0.36/rev 

Chorcferise bending 

5.35 0.70 

4.0 0.62 

Torsion 

9.95 1.30 

9.2 1.48 



TABLE 2.- GAIN MATRIX FOR BELL TABLE 3.- GAIN MATRIX FOR BELL 

PROP ROTOR Gain * 10, Q - 5*.l, PROPROTOR Gain * 10, Q = 6*1, 

4*1., V/fiR = .7. 3*0, V/C2R = .7. 


Controls Controls 












TABLE 4.- GAIN MATRIX FOR BELL 
PROP ROTOR Gain = 10, Q = 6*0, 
3*1, V/fiR = .7. 


TABLE 5.- GAIN MATRIX FOR BELL 
PROP ROTOR W/FLAP CONTROL Gain = 10, 
Q= 5* . 1 , 4*1., V/8R = .7. 


State 

variables 


ontrol 


Controls 


State 

variables 











TABLE 6.- GAIN MATRIX FOR BELL 

PROP ROTOR Gain - 10, Q - 5*.l, 

4*1. , V. = 100. 
knots 


TABLE 7 - GAIN MATRIX FOR BELL 

PROPROTOR Gain - 10, Q * 5*.l, 

4*1. , V, ^ = 250. 

knots 










, .1 , 






r 


'1 


TABLE 8.- 
PRDP ROTOR 

4*1 


GAIN MATRIX FOR BELL 
Gain = 10, Q = 5*.l, 


> V, 


400. 


knots 

Controls 


TABLE 9.- GAIN MATRIX FOR BOEING 
PROPROTOR Gain - 10, Q - 5*. 2, 
10, 1, 10, 1 V/£2R = .7. 

Controls 


State 

variables 

9 lc 

G Is 

e 0 

State 

variables 

e lc 

® Is 

6 0 

3lc 

2.59 

-.29 

-.00 

01c 

-.80 

-.18 

.01 

&lS 

• 45 

1.88 

.03 

Sis 

.12 

-1.10 

- . u5 

he 

.40 

-.02 

.00 

1 he 

1 

4.25 

.01 

-.01 

hs 

-.21 

.48 

-.01 

. 

?ls 

.01 

| 

4.34 

.03 

B 0 

-.06 

.01 

1.38 

So 

.19 

.02 

5.12 


-2.83 

10.05 

.77 

• 

q l 

-8.83 

4.22 

-6.55 

*2 

-2.55 

-5.76 

-4.11 

*2 

16.61 

18.40 

-19.23 

P 

-6.39 

4.39 

-.41 

P 

-4.45 

4.71 

.57 

h 

.06 

-.01 

2.87 

Co 

-.23 

-.03 

4.20 

Sic 

-.02 

O 

. 

c* 

1 

-.00 

Sic 

4.00 

1.07 

.04 

Sis 

2.32 

.22 

-.07 

Sis 

-.79 

4.00 

.04 

Clc 

3.19 

.14 

-.07 

Clc 

2.56 

-4.40 

-.03 

Cls 

.38 

3.01 

-.03 

Cls 

4.07 

2.88 

.13 

So 

.09 

-.04 

.71 

So 

-.24 

-.16 

30.68 

q l 

-3.66 

2.99 

.90 

*1 

-4.07 

7.39 

-1.26 

q 2 

-6.60 

1 

. 

GO 

-4.79 

q 2 

12.75 

9.12 

-4.28 

P 

-4.89 

| 6.93 

-.26 

P 

-2.96 

2.51 

-.24 
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TABLE 12.- GAIN MATRIX FOR BOEING 
PROPROTOR W/FLAP CONTROL Gain = 10, 
Q = 5*. 2, 10, 1, 10, 1, V/flR = .7. 

Controls 


TABLE 13.- GAIN MATRIX FOR 
PROPROTOR Gain = 10, Q = 5*. 

10 ' 1 \nots ' l00 ' 

Controls 


BOEING 
2 , 10 , 


State 

variables 

®ls 

®ls 

0 o 

6 f 

State 

variables 

B 1 c 

^ls 

6 0 

eic 

-.97 

-.06 

.00 

.05 

8 ic 

.76 

-.19 

.00 

3 is 

-.05 

-.88 

-.06 

.08 

61 s 

-.06 

.69 

-.04 

Sic 

4.17 

.09 

.01 

.00 

SlC 

5.03 

.10 

-.09 

Sis 

.17 

4.09 

.03 

-.15 

Sis 

-.06 

5.15 

.03 

60 

.22 

.00 

5.12 

-.09 

BO 

-.02 

-.03 

7.59 

# 

q l 

-1.83 

.62 

-5.56 

32.45 

q l 

.56 

15.31 

-6.76 

*2 

14.22 

20.00 

-18.97 

5.98 

• 

q 2 

25.45 

i 

17.56 

-23.95 

| 

P 

.28 

1.19 

.48 

-1. 39 

• 

P 

-3.74 

5.96 

.47 

So 

-.25 

-.02 

4.20 

.20 

SO 

-.08 

| 

-.00 

3.96 

6 lc 

3.93 

1.01 

.06 

-.17 

Blc 

4.37 

-.86 1 

-.01 

Bis 

-.84 

4.04 

.01 

-.05 

I 

Bis 

.53 

4.12 

-.02 

5 lc 

2.69 

-4.19 

-.01 

.88 

Sic 

.82 

-5.08 

-.10 

Sis 

4.18 

2.72 

.16 

-.07 

Sis 

4.78 

.96 

-.00 

Bo 

-.16 

-.22 

30.66 

-.94 

Bo 

-.31 

-.20 

29.80 

*1 

-1.96 

1.56 

-.75 

3.04 

q l 

-4.72 

10.15 

-.52 

q 2 

12.01 

9.84 

-4.25 

-.03 

q 2 

4.73 

4.36 

-2.03 

P 

-2.39 

.63 

-.29 

-5.89 

P 

-1.02 

-.94 

-.05 


31 



\ , . , j 

4 



1 . v / | 

Jl 


r , S I . < * ■ - s' 

*■ „ | ^ , . . 





















































1 



(a) Rms rotor flap and lag response (closed loop/open loop). 

Figure 7.- Optimal control for Bell rotor at V/OR » 0.7, gain sweep for 

! Q ■ 5*.l, 4*1. 
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Figure 7.- Continued. 
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(d) R«s decoloration for .in, v.rtic.1 bondln, «,) (.!««* 1«P/1>« “>°>» 

Figure 7.- Continued. 


rms occelerotiofl 


! 


-• 4 




0 10 20 30 40 

G«»i 


(f) Ras response of wing and rotor notion (closed loop/open loop). 

Figure 7.- Continued. 





(h) Damping ratio of wing modss 
I-igurfe 7.— Concluded. 
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Figure 8.— Root locus for Bell rotor, gain sweep for Q = 5*.l, 4*1. 
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(a) Rms acceleration of wing and rotor motion (closed loop/open loop). 
Figure 10.— Control of Bell rotor with constraint on rotor motion only, 


Q * 6*1, 3*0. 
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(d) Root locus. 
Figure 10.— Concluded. 
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(c) Damping ratio of wing modes. 
Figure 11.— Continued. 
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(a) Rms acceleration of wing and rotor motion (closed loop/open loop). 

Figure 12.— Optimal control of Bell rotor including wing flap control input, 

Q .= 5* . 1 , 4*1. 
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(c) Damping ratio of wing mr 
Figure 12.— Continued. 
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(a) Rms acceleration, optimum controller (closed loop/open loop). 

Figure 13.— Control of Bell rotor for velocity sweep, with optimum controller 
for each speed and with controller designed for 250 knots (Q = 5*.l, 4*1.). 
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(b) Rms acceleration. 


250 knot controller (closed loop/open loop) 
Figure 13 .- Continued. 
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(c) Rms response, optimum controller (closed loop/open loop). 


















(h) Root locus, 250 knot controller. 
Figure 13.— Continued. 
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(i) Rms response of flap and lag notion. 
Figure 13.— Continued. 
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(j) Rms acceleration of wing vertic 1 bending motion. 

Figure 13.- Concluded. 
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(a) Rms rotor flap and lag response (closed loop/open loop). 

Figure 14.— Optimal control for Boeing rotor at V/OR « 0.7, gain sweep for 

Q = 5*. 2, 10, 1, 10, 1. 
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(e) Rms acceleration of wing ajid rotor motion (closed loop/open loop) . 

Figure 14.— Continued. 
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Cf) Bns response of wing end rotor notion (closed 

Figure 14.^ Continued# 


loop/open loop). 





(g) Steady state response of wing and rotor (closed loop/open loop) . 

Figure 14.- Continued. 
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(h) Damping ratio of wing modes. 

Figure 14.— Concluded. 


1 

I 

























(a) Rms acceleration of wing and rotor motion (closed loop/open loop). 

Figure 18.— Control of Boeing rotor with constraint on wing motion only, 

Q ■ 6*0, 3*1. 
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| (c) Damping ratio of wing modes, 

j Figure 18.— Continued. 
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(c) Damping ratio of wing modes. 
Figure 19.— Continued. 
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(d)'Root locus. 
Figure 19.— Concluded. 
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j (a) Rats acceleration, optimum controller (closed loop/open loop). 


Figure 20.— Control of Boeing rotor for velocity sweep, with optimum controller 
for each speed and with controller designed for 250 knots 
(Q • 5*. 2, 10, 1, 10, 1). 
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(b) Rms acceleration, 250 knot controller (closed loop/open loop). 


Figure 20.- Continued. 
108 




1 


.25 


i 

! 

I 


.20 


.15 



< 


.05 


100 


200 3 < 

V, Knot* 


(e) Damping ratio, optimum < 
Figure 20.- Continue 
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